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Abstract. A new approach to more accurately monitor and 
evaluate transboundary particulate matter (PM) pollution is 
introduced based on aerosol optical products from Korea’s 
Geostationary Ocean Color Imager (GOCI). The area stud- 
ied is Northeast Asia (including eastern parts of China, the 
Korean peninsula and Japan), where GOCI has been mon- 
itoring since June 2010. The hourly multi-spectral aerosol 
optical data that were retrieved from GOCI sensor onboard 
geostationary satellite COMS (Communication, Ocean, and 
Meteorology Satellite) through the Yonsei aerosol retrieval 
algorithm were first presented and used in this study. The 
GOCI-retrieved aerosol optical data are integrated with esti- 
mated aerosol distributions from US EPA Models-3/CMAQ 
(Community Multi-scale Air Quality) v4.5.1 model simu- 
lations via data assimilation technique, thereby making the 
aerosol data spatially continuous and available even for cloud 
contamination cells. The assimilated aerosol optical data are 
utilized to provide quantitative estimates of transboundary 
PM pollution from China to the Korean peninsula and Japan. 
For the period of 1 April to 3 1 May, 2011 this analysis yields 
estimates that AOD as a proxy for PM2.5 or PM10 during 
long-range transport events increased by 117-265% com- 


pared to background average AOD (aerosol optical depth) 
at the four AERONET sites in Korea, and average AOD in- 
creases of 121 % were found when averaged over the entire 
Korean peninsula. This paper demonstrates that the use of 
multi-spectral AOD retrievals from geostationary satellites 
can improve estimates of transboundary PM pollution. Such 
data will become more widely available later this decade 
when new sensors such as the GEMS (Geostationary Envi- 
ronment Monitoring Spectrometer) and GOCI-2 are sched- 
uled to be launched. 


1 Introduction 

Transboundary air pollution is an important issue in Asia 

(Arndt et al., 1998; Holloway et al., 2002; Yamaji et al., 

2012). In particular, the long-range transport (LRT) of air 

pollutants from China (source) to the Korean peninsula and 

Japan (receptors) due to persistent westerlies over North- 

east Asia has been a central issue in this region (Arndt et 

al., 1998; Song and Carmichael, 2001; Yamaji et al., 2012). 

Transboundary air pollution is a serious issue not only in 
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Asia, but also in Europe and North America, causing a num- 
ber of different international problems such as acid precipi- 
tation, rising levels of hemispheric ground-level ozone, and 
regional visibility impairment (Levy II and Moxim, 1987; 
Brankov et al., 2003; Holloway et ah, 2003; Tulet et ah, 
2003; Park et al., 2004; Wang et al., 2009). To manage trans- 
boundary air pollution between countries, the Convention on 
Long-Range Transboundary Air Pollution (CLRTAP) under 
the UNECE (United Nations Economic Commission for Eu- 
rope) was launched to lessen air pollution through the LRT 
in 1979 (http://www.unece.org/env/lrtap/lrtap_hl.html). The 
EMEP (European Monitoring and Evaluation Programme) is 
a scientific effort established to support the CLRTAP. 

Quantitative evaluation of transboundary air pollution re- 
mains a difficult task and currently cannot be done based 
on observations alone due to a lack of explicit source sig- 
natures and lack of an appropriate observation system (Mal- 
colm et al., 2000; Brankov et al., 2003; Tulet et al., 2003; 
Park et al., 2004; Schaub et al., 2005; Wang et al., 2009). 
Most current evaluations of transboundary transport of pol- 
lutants are based on chemistry-transport model (CTM) sim- 
ulations (Levy II and Moxim, 1987; Arndt et al., 1998; Hol- 
loway et al., 2002; Wang et al., 2009; Yamaji et al., 2012). By 
controlling the emissions from a source country (e.g., China) 
in the CTM simulations, transboundary air pollution from the 
source country to receptor regions has been evaluated. How- 
ever, these model-based evaluations have significant limita- 
tions because of the large uncertainties in the emission rates 
of the air pollutants from the source country (or region), 
the highly nonlinear nature of atmospheric chemistry, and 
uncertainties related to the CTM and meteorological mod- 
els (MM). For example, both atmospheric secondary parti- 
cle formation and primary particle (e.g., mineral dust) gen- 
eration/transport are difficult to correctly simulate (Song and 
Carmichael, 2001; Song et al., 2005; Volkamer et al., 2006; 
Robinson et al., 2007). These uncertainties are a main reason 
for disagreement in source and receptor estimates made by 
different countries, which has led to distrust of the outcomes 
of model-based assessments. 

Current efforts to reduce the uncertainty in estimates of 
long-range transport of pollutants by integrating models and 
observations through data assimilation are limited by our cur- 
rent observational systems. Surface-based observations are 
too sparse to constrain the system, while current satellites are 
unable to retrieve aerosol concentrations. Instead, satellite re- 
trievals of AOD are used to help constrain modeled AOD and 
the modeled relationship between AOD and surface-level PM 
(particulate matter) can then be used to estimate surface-level 
PM concentrations (Wang and Christopher, 2003; Liu et al., 
2005; van Donkelaar et al., 2006; Schaap et al., 2009; Tsai et 
al., 201 1). This approach is currently limited by the fact that 
AOD from polar-orbiting satellites produce retrievals only 
once per day (Kaufman et al., 1997; Diner et al., 2001) and 
current geostationary satellites lack “multi-spectral” aerosol 
data (Knapp et al., 2002; Wang et al., 2003). For example, 
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the USA GOES series (GOES 8, 13), the Japanese GMS and 
MTSAT, and European SEVIRI have been meteorological 
sensors with only one to two visible channels. The history 
of GEO satellite sensors for the retrieval of the aerosol opti- 
cal properties is summarized in Table 1 . 

This situation is changing as new geostationary satel- 
lite based sensors are becoming available and will be ex- 
panded by upcoming planned missions. The first hourly 
multi-spectral aerosol data observed from a geostationary 
satellite sensor, GOCI (Geostationary Ocean Color Imager) 
is used in this paper to demonstrate how this new informa- 
tion can be used to improve our estimates of LRT of pol- 
lutants. Specifically, the level 1 radiance data from GOCI, 
onboard COMS (Communication, Ocean, and Meteorology 
Satellite), is converted into level 2 AOD products via the 
“Yonsei aerosol retrieval algorithm” (the details of the al- 
gorithm will be discussed in the Methods section) (Lee et 
al., 2010). COMS/GOCI was launched in June 2010, and the 
COMS/GOCI data center began to release radiance data in 
April 2011, after an initial operation test (IOT) period. This 
satellite sensor was designed for monitoring ocean colors 
around the Korean peninsula, but it can also provide AOD 
data for a strategically important region, i.e., Northeast Asia. 
The domain covered by the GOCI sensor shown in Fig. 1 
is 2500 km x 2500 km and includes eastern parts of China, 
the Korean peninsula, and Japan. The GOCI-retrieved AOD 
products were integrated with CMAQ-simulated AOD data 
to produce the more accurate aerosol data set used to investi- 
gate the impacts of transboundary PM pollution between the 
source and receptor countries in Northeast Asia. 

2 Methods 

The procedure to produce GOCI-constrained AOD product 
combines the GOCI level 2 product with AOD predictions 
from the US EPA Model-3/CMAQ (Community Multi-scale 
Air Quality) v4.5.1 model simulations. The CMAQ model 
simulations were performed over East Asia, with emission 
data prepared appropriately for the model domain and mete- 
orological fields produced from the PSU/NCAR MM5 (the 
fifth generation meteorological model) simulations. The fi- 
nal AOD product is spatially continuous and available even 
for locations where GOCI retrievals are not available (e.g., 
cloud-contaminated points, sun-glint areas, desert, and snow- 
covered land). The use of GOCI-retrievals overcomes the 
“temporal limitation” typical in analyses with LEO (low 
Earth orbit) aerosol data and the assimilation overcomes the 
“spatial limitation” that has been typical in analyses based 
only on satellite data. The assimilated AOD data set using the 
GOCI data is utilized to monitor and evaluate transboundary 
particulate pollution from China to the Korean peninsula in 
this study. 
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Table 1 . Brief descriptions of aerosol optical property retrievals from some geostationary satellites. 


Satellite & 
Sensor 

GMS5 1 
VIS SR 7 

MTSAT-1 R 2 
JAMI 8 

MSG 3 

SEVIRI 9 

GOES 4 

GOES-8 

COMS 5 

GOCI 10 

GOES-R 

ABI 11 

GEO-KOMPSAT 6 
GEMS 12 & ABI & 
GOCI-2 

Sensor type 

MET 13 sensor 

MET sensor 

MET sensor 

MET sensor 

Ocean color imager 

MET sensor 

Environ. & ocean 
color imager 

Temporal resolution 

lh 

30 min 

15 min 

15 min 

lh 

TBD 17 

1 h (GEMS) 

Spatial resolution 

a 1.25km/ b 5km 

1.25 km/ 5 km 

1 km/ 3 km 

1 km/ 4 km 

0.5 km/ 0.5 km 

0.5 km /2 km 

0.25 km/ 1km 

Number of bands 

1 VIS 14 / 3 IR 16 

1 VIS/4 IR 

2 VIS/2NIR 15 / 
8 IR 

1 VIS MIR 

6 VIS/2NIR 

2 VIS / 4 NIR / 
10 IR 

13 VIS/NIR 

Covering region 

East Asia 

& western Pacific Ocean 

East Asia 

& western Pacific Ocean 

Central 

Europe 

Eastern 

US 

Northeast 

Asia 

US 

East, Southeast 
& Central Asia 

Launch year 

Mar 1995 

Feb 2005 

Aug 2002 

Apr 1994 

Jun 2010 

Scheduled to be 
launched in 2015 

Scheduled to be 
launched in 2018 

References 

Masuda et al. (2002); 
Wang et al. (2003) 

Kim et al. (2008) 

Popp et al. (2007) 

Christopher et al. (2002); 
Knapp et al. (2002, 2005) 

Lee et al. (2010); 
this study 

Laszlo et al. (2008) 



1 GMS: Geostationary Meteorological Satellite, 2 MTSAT-1R: Multi-functional Transport Satellite, 3 MSG: Meteosat Second Generation, 4 GOES: Geostationary Operational Environmental Satellite, 5 COMS: Communication, Ocean, and Meteorological 
Satellite, 6 GEO-KOMPSAT: Geostationary Earth Orbit KOrea Multi-Purpose SATellite, 7 VISSR: Visible and InfraRed Spin Scan Radiometer, 8 JAMI: Japanese Advanced Meteorological Imager, 9 SEVIRI: Spinning Enhanced Visible and InfraRed Imager, 
10 GOCI: Geostationary Ocean Color Imager, 1 1 ABI: Advanced Baseline Imager, 12 GEMS: Geostationary Environment Monitoring Spectrometer, 13 MET: Meteorological, 14 VIS: Visible, 15 NIR: Near InfraRed, 16 IR: InfraRed, l7 TBD: To Be 
Determined, a Resolution for the visible channels, ^Resolution for the infrared channels. 


2.1 Meteorology and chemistry-transport model 
simulations 

The PSU/NCAR MM5 and US EPA Model-3/CMAQ model 
v4.5.1 were respectively used as the meteorological model 
(MM) and chemistry-transport model (CTM) in this study 
(Byun and Schere, 2006; Park et ah, 2011). The accuracy 
of the MM5 simulations was enhanced by nudging with 
the observed meteorological data (Grell et al., 1995). Ini- 
tial conditions (ICs) and boundary conditions (BCs) for the 
MM5 simulations were updated continuously with the Na- 
tional Centers for Environmental Prediction (NCEP) reanal- 
ysis data. The latest remotely sensed surface wind data re- 
trieved from the Advanced SCATterometer (AS CAT), which 
was launched on 19 October 2006 onboard the Meteoro- 
logical Operational satellite A (MetOp-A) was used for 
the four-dimensional data assimilation (FDDA), since dur- 
ing regional-scale, long-range transport events the signif- 
icant fractions of air pollutants are frequently transported 
through the marine boundary layer (MBL) between countries 
in Northeast Asia (e.g., through the MBL over the Yellow Sea 
between China and Korea and through the MBL over the East 
China Sea between China and Japan; regarding the locations 
of the Yellow Sea and the East China Sea, refer to Fig. 1) 
(Daum et al., 1996; Rajeev et al., 2000; Song et al., 2007). 

The schemes chosen for CMAQ model simulations were 
the same as those used in previous research (Park et al., 
2011). The SAPRC-99 chemical mechanism and AER04 
aerosol module were selected for simulating gas-phase 
chemistry and atmospheric aerosol thermodynamic/dynamic 
processes, respectively. More details of the MM5 and CMAQ 
model simulations were described in a previous publication 
(Park et al., 2011). 

The CMAQ model simulations were carried out using the 
best available emission inventories over East Asia that have 
been made more sophisticated with, for example, the appli- 
cations of monthly variations of anthropogenic emissions, 
the uses of an operational dust forecast model at the KMA 
(Korea Meteorological Administration) and the NCAR (Na- 



Fig. 1 . Domains of CMAQ model simulations (in blue) and 
COMS/GOCI sensor (in red). Also shown are five AERONET sites 
in the receptor areas (Korea and Japan). For the sake of analytical 
convenience, three subregions (A), (B), and (C) in the Korean penin- 
sula are defined; these represent the northern, central, and southern 
parts of the Korean peninsula, respectively. 


tional Center for Atmospheric Research) FINN (Fire Inven- 
tory from NCAR) biomass-burning (BB) emissions (Wied- 
inmyer et al., 2011). The efforts are summarized in Table 2. 
For the consideration of anthropogenic pollutants emissions, 
such as NO x , SO 2 , CO, black carbon (BC), NMVOCs, etc., 
Intercontinental Chemical Transport Experiment-Phase B 
(INTEX-B), Regional Emission inventory in Asia (REAS), 
and Clean Air Policy Support System (CAPSS) were applied 
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as the best suitable emission data over China, Japan, and 
the Korean peninsula, respectively. These emission invento- 
ries originally adjusted to the year 2006 were updated for 
the year 2011 modeling period, using the projection fac- 
tors based on the Greenhouse gas and Air pollution INterac- 
tions and Synergies (GAINS) modeling emissions scenarios 
(http://gains.iiasa.ac.at/index.php/gains-asia). In the mean- 
while, in order to consider isoprene and monoterpene emis- 
sions over the modeling domain, the isoprene emissions from 
Model of Emissions of Gases and Aerosols from Nature 
(MEGAN) and MOdel for Hydrocarbon emissions by the 
CANopy (MOHYCAN) were utilized for the CMAQ model- 
ing input data (Muller et al., 2008). Monoterpene emissions 
were also estimated from the isoprene emissions, using the 
ratios of monoterpene to isoprene from Global Emissions In- 
ventory Activity (GEIA). 

In particular, this modeling period (a spring season in 
Northeast Asia) could not be explained without considering 
the effects of BB and mineral dust emissions. FINNvl is a 
BB emission inventory annually provided from the year 2002 
by NCAR (http://bai.acd.ucar.edu/Data/fire/). This FINNvl 
has produced the BB emission data with a high horizontal 
resolution of ~ 1 km x 1 km over the entire globe, being ap- 
plied for variable regions from local to global scales (Wied- 
inmyer et al., 201 1). In this study, the FINNvl was used for 
the BB emissions. The CMAQ v4.5.1 model does not in- 
clude any dust modeling modules related to occurrence and 
transport of mineral dust events. However, without taking the 
influence of mineral dust events into account, accurate esti- 
mation of the aerosols effects observed by long-range trans- 
ports in the Korean peninsula and Japan would not be possi- 
ble. Therefore, the outbreak and transport of Asian mineral 
dust were considered by Asian Dust Aerosol Model (ADAM) 
simulations. The ADAM is currently operational model for 
forecasting Asian mineral dust events at the KMA and the 
details of the ADAM simulations were reported previously 
by e.g., Park and Lee (2004), Kim et al. (201 1), and Park et 
al. (2011). 

The domain of the CMAQ model simulation is also shown 
in Fig. 1. The CMAQ domain (blue color) covered the re- 
gion of approximately 95-145° E, 20-50° N, which includes 
COMS/GOCI domain (red color) composed of Northeast 
China, the Korean peninsula, and Japan. In this study, the 
CMAQ domain was chosen to be larger than the GOCI 
domain, since the air quality inside the GOCI domain is 
strongly influenced by the influxes from the outer areas with 
persistent westerlies in East Asia. The horizontal grid spac- 
ing of the CMAQ modeling was 30 km x 30 km and the 
vertical domain was distributed in 14 layers following o- 
coordinate from lOOOhPa to 70hPa. The MM5-CMAQ sim- 
ulation period started at 00:00 UTC, 1 April to 23:00 UTC, 
31 May 2011 (with a week-long CMAQ spin-up period), 
since the first batch of the level 1 radiance data from the 
COMS/GOCI data center was released for this period. 
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2.2 Remote-sensing data 
2.2.1 COMS/GOCI data 

Sensors onboard LEO satellites, such as MODIS (Moderate 
Resolution Imaging Spectroradiometer) and MISR (Multi- 
angle Imaging Spectroradiometer), have ordinarily been uti- 
lized to obtain information on aerosol optical properties 
(AOPs), but their monitoring has taken place at most once 
per day. To more accurately monitor the transboundary 
LRT events of air pollutants, hourly (or highly temporally 
resolved) monitoring data from geostationary Earth orbit 
(GEO) sensors are needed instead of the temporally lim- 
ited LEO measured data. As mentioned previously and also 
presented in Table 1, MSG, GOES, and MTSAT-1R sen- 
sors have been partly monitoring the AOPs over Europe, 
US, and Japan, respectively (Knapp et al., 2005; Popp et 
al., 2007; Kim et al., 2008), but these meteorological sen- 
sors are equipped with only one visible channel, which pos- 
sibly makes the quality of the AOPs from these meteoro- 
logical sensors less accurate. In contrast, GOCI sensor on- 
board geostationary COMS has six visible bands at 412, 
443, 490, 555, 660, 680 nm, which can make advances by 
more accurately calculating the AOPs, compared to other 
meteorological GEO satellite sensors. In addition, atmo- 
spheric compensations, such as cloud-screening make the 
numerical accuracy of the AOPs more improved, because 
the COMS/GOCI has 2 more bands at 745 nm and 865 nm 
(i.e., near IR). The GOCI sensor covers Northeast China, the 
Korean peninsula, and Japan, and the products can be gen- 
erated with 500 m x 500 m spatial resolution. GOCI level 1 
data during daytime can be obtained eight times from 09:00 
to 16:00 LST (local standard time) (refer to http://kosc.kordi. 
re.kr/oceansatellite/coms-goci/introduction.kosc). 

The hourly GOCI level 1 data is released somewhat slowly 
after the IOT and data calibration procedures. The GOCI 
level 1 radiance is then converted into level 2 AOD data 
using a multi-spectral algorithm being developed by Jhoon 
Kim and coworkers at Yonsei University, Korea (Lee et al., 
2010). The level 2 AOD data retrieved with a pixel size of 
6 km x 6 km was used in this study. The retrieval algorithm 
was originally based on the NASA MODIS (i.e., LEO) al- 
gorithm, but it has evolved to date to include the aerosol 
models optimized for East Asia and to consider the geometry 
from GEO satellite instruments such as GOCI, GOCI-2, and 
GEMS (Geostationary Environment Monitoring Spectrome- 
ter) (both the GOCI-2 and GEMS sensors are scheduled to 
be launched in 2018). The Yonsei aerosol retrieval algorithm 
includes the capability to classify aerosol types based on ab- 
sorption tests and size information (Kim et al., 2007; Lee et 
al., 2010). The evaluation of the current AOD data from the 
algorithm shows that the slope of the linear regression line 
is close to 1.03 over ocean, but 1.15 over land when com- 
pared to values from AERONET. Although the detailed ex- 
planation and validation procedures of this Yonsei aerosol 
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Table 2. Emission inventories used in the CMAQ model simulations. 


Region 

Emission inventory 

Species 

China and North Korea 

INTEX-B 1 

S0 2 , NO x , NMVOCs, CO, C0 2 , BC, OC 


EDGAR 2 

nh 3 

South Korea 

CAPSS 3 

S0 2 , NO x , NMVOCs, CO, NH 3 


INTEX-B 

C0 2 , BC, OC, Chemical split factors of NMVOCs 

Japan 

REAS 4 

S0 2 , NO x , NMVOCs, CO, C0 2 , NH 3 , BC, OC 

Entire domain 

GAINS 5 emission scenarios 

a Interannual variations of CO, NO x , VOC, S0 2 , NH 3 , C0 2 , CH 4 , PM 2 5 


FINN 6 fire emissions 

CO, NO, N0 2 , S0 2 , NH 3 , NMVOCs, BC, OC, PM 2 5 , PM 10 


MEGAN 7 +MOHYC AN 8 

Isoprene 


GEIA 9 

^Ratios of monoterpene to isoprene 


1 INTEX-B: Intercontinental Chemical Transport Experiment-Phase B, 2 EDGAR: Emission Database for Global Atmospheric Research, 3 CAPSS: Clean Air Policy Support 
System, 4 REAS: Regional Emission inventory in ASia, 5 GAINS: Greenhouse gas and Air pollution INteractions and Synergies, 6 FINN: Fire Inventory from NCAR, 

7 MEGAN: Model of Emissions of Gases and Aerosols from Nature, 8 MOHYCAN: MOdel for Hydrocarbon emissions by the CANopy, 9 GEIA: Global Emissions Inventory 
Activity, a Emission increase (or decrease) rates from 2006 to 201 1 for the air pollutants were obtained from the GAINS emission scenarios, b Monoterpene emissions are not 
included in MEGAN+MOHYCAN emissions. Therefore, the ratios of monoterpene to isoprene in the biogenic emissions were obtained from the GEIA emission inventory. 


retrieval algorithm were given in previous publications (Lee 
et al., 2010, 2012), a brief introduction on the algorithm is 
made here. The retrievals of the AOPs were performed with 
top of atmosphere spectral reflectances. To obtain the AOPs 
over the regions including clouds, clouds are screened by 
spatial variability and threshold tests at each pixel. Even so, 
there are limitations to the complete removal of clouds and 
the perfect retrieval of surface reflectance due to the lack of 
IR bands on GOCI. Therefore, the Yonsei aerosol retrieval al- 
gorithm eliminates clouds utilizing a meteorological imager 
(MI) that has 4 IR channels. The MI sensor is also onboard 
COMS with GOCI. In addition, the Yonsei aerosol retrieval 
algorithm is composed of two different retrieval logics that 
are applied to “continent” including near-coastal turbid wa- 
ter areas and “ocean”, respectively. 

Figure 2 presents one example: a comparison of the 
Terra/MODIS-retrieved and Aqua/MODIS-retrieved AODs 
with GOCI-retrieved AODs at approximately 10:00 LST and 
13:00 LST on 10 April and 12 April, respectively. As shown 
in Fig. 2, there are agreements in the spatial AOD distri- 
butions from the two MODIS and GOCI sensors. As dis- 
cussed previously, at most two snapshots of information per 
day from the LEO sensors (i.e., MODIS) can be obtained, 
whereas eight snapshots per day from the GEO sensor (i.e., 
GOCI) are available. This is an obvious “temporal” advan- 
tage that can be achieved from the GEO sensor. Furthermore, 
as shown in Fig. 2, the AOD data from the GOCI sensor 
shows a better spatial coverage, compared to those from the 
MODIS (see white bands along the north-south swath in the 
MODIS data). This lack of the MODIS data was caused by 
sun-glint over the ocean areas (details about the sun-glint ef- 
fects are discussed further in Appendix A). By contrast, such 
lack of the data is not found in the GOCI data. This is an- 
other “spatial” advantage that can be achieved with the GEO 
sensor. 


(a) MODIS-retrieved AOD (b) GOCI-retrieved AOD 


120E 130E 140E 120E 130E 140E 



Fig. 2. Snapshots of (a) Terra and Aqua/MODIS-retrieved AODs 
and (b) GOCI-retrieved AOD over Northeast Asia at 10:00 LST 
(Terra/MODIS) and 13:00 LST (Aqua/MODIS) on 10 April and 
12 April. 
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2.2.2 AERONET 

Many studies have used the ground-based NASA AERONET 
data to verify modeling results and satellite monitoring data 
(e.g., Chin et al., 2002, 2004; Yumimoto and Takemura, 
2011). The AERONET data measured by sun photometers 
were utilized in this study to investigate the accuracy of 
CMAQ model simulations and GOCI AOP data. 

In this manuscript, the data obtained from five AERONET 
observation sites (i.e., Baengnyeong, Yonsei Univer- 
sity (YonseiUniv), GIST (GwangjuGIST), Gosan 
(Gosan SNU), and Osaka) in the Korean peninsula and 
Japan were selected in order to analyze long-range transport 
effects over Northeast Asia. The AERONET AOD level 2.0 
data used for comparison was obtained from the official 
NASA web site at http://aeronet.gsfc.nasa.gov. The only 
exception was the AERONET AOD data for the Yonsei 
University (Seoul) site, which were not calibrated. Before 
use, they were upgraded by eliminating some sporadic cloud 
contamination of the initial high AOD data (AOD >1.0) 
using the meteorological data from the automatic weather 
stations (AWS) colocated with Cimel sun photometers. 
The uncertainty and possible anomalies in the Cimel sun 
photometer AOD data at the Yonsei University site were also 
checked by the Laboratory for Terrestrial Physics at NASA 
Goddard Space Flight Center. AOD uncertainties of 0.01 to 
0.04 were estimated, mainly due to calibration. 

2.3 AOD calculations from CMAQ model simulations 

Concentrations of nitrate, sulfate, ammonium, organic 
aerosols (OAs), BC, sea-salt aerosols, and mineral dust are 
provided as results of the CMAQ model simulations. This 
aerosol composition can be converted into the AOPs includ- 
ing AOD by considering the mass extinction efficiencies 
(MEEs) depending on the scattering and absorbing charac- 
teristics for each aerosol composition and the hygroscopicity. 
First, the absorption, scattering, and extinction coefficients 
O^abs? ^scat? and a ex t) were calculated by the following for- 
mulas: 

^ext (^) — ^abs^) T - ° r scat(^)? (1) 


cr abs (a. : Mm- 1 ) = {(1 ~OHj W) • fty W • • [C]y}, (2) 

f^scat | r - ■ Mm ) = et; {wjM-Pijfr)-fij($W'[C]ij}, (3) 

where i and j represent the particulate species i at the size 
bin or particle mode j ; fij (RH) the hygroscopic enhance- 
ment factor as a function of relative humidity (RH); [Clo- 
the concentrations of particulate species such as (NHz^SCU, 
NH 4 NO 3 , BC, OAs, mineral dust, and sea-salt aerosols; and 
c Oij (A.) and fyj (X) the single scattering albedo and MEE of 
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particulate species i at the size bin or particle mode j, re- 
spectively. The terms related to the MEEs (/3; ; ) and the hy- 
groscopicity (//j(RH)) have been defined differently in sev- 
eral studies, e.g., d’Almeida et al. (1991), Malm et al. (1994), 
Chin et al. (2002), Malm and Hand (2007), and Park et 
al. (201 1). According to the previous investigations, the way 
to compute the extinction coefficients ( a ext ) has been mainly 
classified into two categories. The first method is to use em- 
pirical equations (e.g., Malm et al., 1994; Malm and Hand, 
2007; Song et al., 2008; Park et al., 2011). The other method 
is to use the equations obtained based on the Mie theory (e.g., 
d’Almeida et al., 1991; Chin et al., 2002). In this study, the 
methods from Chin et al. (2002) and Park et al. (201 1) were 
mainly utilized and compared. 

Based on the particulate composition and meteorolog- 
ical parameters obtained from the MM5-CMAQ model 
simulations, the AOD was calculated at a wavelength of 
550 nm (AOD 550 ), using the conversion algorithms of Chin 
et al. (2002), Park et al. (201 1), and others (d’Almeida et al., 
1991; Malm et al., 1994; Hess et al., 1998; Malm and Hand, 
2007): 

AOD 550 = [ (Text (550 nm) dz , (4) 

Jo 

where z represents the vertical height. For example, the 
CMAQ-simulated AOD products calculated from Park et 
al. ( 2011 ) usually showed underestimations when com- 
pared to the AERONET AOD data. However, when Chin 
et al.’s (2002) conversion algorithm was applied, the AOD 
values were usually larger than those calculated by Park et 
al.’s (2011) method. These results will be discussed further 
in a later section. 

2.4 Data assimilation 

In several previous studies, optimal interpolation (OI) tech- 
nique with a Kalman gain matrix has been employed to im- 
prove the predictions of AOD and aerosol mass distributions 
with the combination of the modeling results and monitor- 
ing data (Collins et al., 2001; Yu et al., 2003; Adhikary et 
al., 2008; Chung et al., 2010; Park et al., 2011). For exam- 
ple, Adhikary et al. (2008) showed some enhancements in 
the estimations of the AOD values over Asia by combin- 
ing the STEM-modeled and MODIS-retrieved data through 
the OI technique. The accuracy of the AOD values was also 
greatly improved after the application of the OI method using 
CMAQ-simulated and MODIS-retrieved AOD data (AOD m 
and AOD 0 in Eq. 5) (Park et al., 201 1). 

AOD m = AOD m + K[AOD 0 - H(AOD m )], (5) 

K=BH r (HBH r +0)~ 1 , ( 6 ) 

where H denotes a linear operator for interpolation from the 
modeling grid to the satellite retrieval grid. B and O represent 
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the error covariances of the modeling and satellite (observa- 
tion) fields, respectively. 

In this study, the 01 technique was again utilized for the 
assimilation of the AOD data from the GOCI sensor into the 
CMAQ model. The Kalman gain matrix (K in Eq. 6) was 
obtained by the covariance of modeled and observed fields 
(B and O in Eq. 6) which is a function of free parameters. 
The yearly fixed free parameters were mainly used in sev- 
eral previous studies (Collins et al., 2001; Yu et al., 2003; 
Adhikary et al., 2008; Chung et al., 2010), whereas Park et 
al. (2011) proposed applying seasonal variations of free pa- 
rameters. The best free parameters for four seasons suggested 
in Park et al. (2011) were used here. Detailed explanations 
of the OI technique including the selection of the free pa- 
rameters were made in Park et al. (2011). The assimilated 
AOD data produced by combining the CMAQ-simulated and 
the GOCI-retrieved AODs were used to investigate the long- 
range PM transport events in this study. 

2.5 Long-range transport of particulate matter 

Transboundary transports of the AOD plumes from the China 
continent to the Korean peninsula and Japan were the fo- 
cus of this study. However, defining the LRT events over 
Northeast Asia is not easy, mainly due to complexities in 
air mass movements. We considered large-scale AOD plumes 
that caused the peak AOD values larger than 0.4 at the recep- 
tor sites (“peak AOD filter”). Such LRT AOD plumes were 
further traced back via NOAA HYSPLIT (HYbrid Single- 
Particle Lagrangian Integrated Trajectory) model analysis 
(Draxler and Hess, 1998) from the points of time when the 
peak AOD values occurred. The HYSPLIT model simula- 
tions were performed with meteorological data, which were 
taken from the GDAS (Global Data Assimilation System) of 
NCEP. When the backward trajectories from the occurring 
points of the peak AOD values had passed through China 
within 48-72 h, the AOD plumes were finally determined to 
be LRT plumes (“backward trajectory filter”). A LRT case 
shown in Lig. 3 was decided by using both of the filters. 
Lor example, 72 h HYSPLIT backward trajectory analysis at 
three AERONET receptor sites in Korea for the LRT case 
in Lig. 3 is presented in Lig. 4. As shown, all the backward 
trajectories were passed through China. By contrast, Lig. 5 
shows an example of a non-LRT case. Both the LRT and non- 
LRT cases will be discussed further in Sect. 3.1. 

In order to investigate the LRT of particulate matter over 
Northeast Asia, the AOD values instead of PM concentra- 
tions were used. In other words, the AOD has been regarded 
as a proxy quantity for PM concentration in this study. In this 
context, relationships between PMio (or PM2.5) and AOD 
can be important. The relationships are affected by various 
parameters, such as height of the PBL (planetary boundary 
layer), RH, vertical distribution of aerosols, seasonality, di- 
urnal variation, geographical information, particle composi- 
tions, and size distributions (Wang and Christopher, 2003; 



Fig. 3. Spatial and temporal distributions of (a) CMAQ-simulated 
AOD (calculated via Chin et al.’s (2002) conversion algorithm), 
(b) GOCI-retrieved AOD (retrieved via the Yonsei aerosol retrieval 
algorithm (Lee et al., 2010)) and (c) assimilated AOD with two- 
hour resolution over Northeast Asia at 09:00, 11:00, 13:00, and 
15:00 LST on 10 April 2011. 


Engel-Cox et al., 2004; Liu et al., 2005; van Donkelaar et al., 
2006; Choi et al., 2009; Schaap et al., 2009; Tsai et al., 201 1). 
Several studies have estimated two-dimensional PM 10 from 
the relationships between ground-based PM10 and satellite- 
retrieved AOD (Wang and Christopher, 2003; Engel-Cox et 
al., 2004; Liu et al., 2005; Schaap et al., 2009; Tsai et al., 
201 1) or the ratios of CTM-simulated PM 10 to AOD (Liu et 
al., 2004; van Donkelaar et al., 2006; Choi et al., 2009). 


3 Results and discussions 

The hourly assimilated AOD data was used to identify and 
analyze the PM LRT events over Northeast Asia. Based on 
this, the variations of the AOD values during LRT events in 
comparison with those during non-LRT events were analyzed 
on a more quantitative basis. 

3.1 Long-range transport of aerosols from China to 
Korea 

A representative snapshot of an LRT event of particulate 
matter from China to the Korean peninsula, Japan, and the 
northwestern Pacific Ocean is presented in Fig. 3. It can 
be seen that a huge AOD plume (AOD >0.6) was gener- 
ated in the greater Beijing area and Shandong province on 
10 April and was then transported over long distances to the 
Korean peninsula following a low pressure trough. GOCI 
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Backward trajectories 


120E 130E 140E 



Fig. 4. NOAA HYSPLIT 72 h backward trajectory analysis at three 
Korean AERONET sites (Baengnyeong, YonseiUniv (Seoul), and 
Gwangju GIST) with 2h resolution at 09:00, 11:00, 13:00, and 
15:00 LST on 10 April 2011. 


retrievals with one-hour intervals from 09:00 to 16:00 LST 
(eight snapshots per day) on 10 April show high aerosol load- 
ings (AOD > 1 .2) over the Shandong peninsula and Shanghai 
regions that extend over Korea. The plumes extend farther 
out over the Yellow Sea and slowly shift to lower latitudes as 
time progresses. In Fig. 3, only two-hour resolution results 
are presented. It is important to note that GOCI AOD val- 
ues are not retrieved in regions where clouds exist. As shown 
in Fig. 3, cloud-contaminated pixels appear white, which in- 
dicates that no data is available at these pixels. The GOCI 
retrievals, when assimilated into the CMAQ model, produce 
complete spatial distributions where AOD is estimated in all 
pixels (see the third column of Fig. 3). Therefore, the assim- 


(a) CMAQ AOD (b) GOCI AOD (c) Assimilated AOD 

120E 130E 140E 120E 130E 140E 120E 130E 140E 



0.0 0.3 0.6 0.9 1.2 1.5 


Fig. 5. As Fig. 3, for a non-LRT event on 12 April 2011. This case 
was determined as non-LRT because the main AOD plumes were 
transported into over northwestern Pacific Ocean, not over the Ko- 
rean peninsula and/or Japan. 


dated AOD fields are much improved in comparison to the 
predicted (a priori) CMAQ AOD distributions (first column 
of Fig. 3). The a priori CMAQ AOD values have a tendency 
to be underestimated compared to the GOCI-retrieved AOD 
values, as shown in Fig. 3. This underestimation is believed 
to be related to uncertainties in the emission rates of gaseous 
and particulate pollutants and the CMAQ model simulations, 
as well as the uncertainties in the assumptions made during 
the AOD estimation from the CMAQ modeling results. 

In this analysis, the AOD calculation via Chin et 
al.’s (2002) conversion method was used. However, sev- 
eral other estimations were also tested (e.g., Malm et al., 
1994; Chin et al., 2004; Park et al., 2011). One example 
is shown in Fig. 6 and Table 3, which were produced by 
Park et al.’s (201 1) conversion method. It was found that the 
AOD values were underestimated because the MEE values 
and hygroscopic enhancement factor (/; 7 (RH)) of Park et 
al.’s (2011) conversion method were reported to be in gen- 
eral smaller than those of Chin et al.’s (2002) method. After 
assimilation both the magnitudes and spatial distributions are 
significantly improved. The AOD LRT event shown in Fig. 3 
is only an example. Seven to fourteen small- and large-scale 
AOD LRT events occurred over the two-month period (from 
1 April to 31 May 2011). 

Furthermore, the non-LRT events were also identified. 
During the non-LRT periods, air masses were generally 
stationary. Although the influx of air pollutants emitted from 
the highly polluted China regions to the Korean peninsula 
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Table 3. Comparison analysis of the assimilated AOD for LRT and non-LRT events at five AERONET sites during the entire modeling period 
of 1 April to 3 1 May 2011. 




Non-LRT events 



LRT events 


Data fraction of 

AOD increase 

Ratio of AOD increase 










LRT events 

during LRT 

during LRT 










Glrt in Eq. 8) 

events 3 

events 4 


N l 


Avg. 

Max. 

N 

o 

Avg. 

Max. 




Baengnyeong 

849 

0.11 

0.16 

0.95 

606 

0.36 

0.59 

2.15 

0.42 

0.43 

2.59 

Yonsei_Univ (Seoul) 

821 

0.15 

0.22 

0.77 

634 

0.29 

0.49 

2.13 

0.44 

0.28 

1.28 

GwangjuGIST 

1024 

0.13 

0.20 

1.10 

431 

0.31 

0.47 

1.89 

0.30 

0.27 

1.36 

GosanSNU 

967 

0.12 

0.17 

1.09 

488 

0.24 

0.41 

1.53 

0.34 

0.24 

1.37 

Osaka 

1247 

0.11 

0.16 

1.16 

208 

0.21 

0.39 

1.23 

0.14 

0.23 

1.41 

Total 


0.12 

0.18 

1.16 


0.31 

0.49 

2.15 

0.33 

0.31 

1.68 


The calculations were conducted with the assimilated AOD products. The AOD data assimilations were carried out with the GOCI-retrieved AOD products (retrieved via the Yonsei algorithm (Lee et 
al., 2010)) and the CMAQ-simulated AOD products (calculated via Park et al.’s conversion algorithm (Park et al., 201 1)). 

1 The number of data for non-LRT or LRT events for the entire modeling period; 

2 standard deviation; 

3 AOD increased by LRT events above the average background AOD calculated during non-LRT periods; 

4 the ratios were estimated by ((AODlrt - AOD non _LRT)/AOD non _LRr). Here, AODlrt and AOD non _LRT indicate the average AOD calculated during the LRT and non-LRT events, respectively. 


Table 4. Comparison analysis of the assimilated AOD for LRT and non-LRT events at five AERONET sites during the entire modeling period 
of 1 April to 3 1 May 2011. 




Non-LRT events 



LRT events 


Data fraction of 

AOD increase 

Ratio of AOD increase 










LRT events 

during LRT 

during LRT 










Glrt in Eq. 8) 

events 3 

events 4 


N l 

O' 2 

Avg. 

Max. 

N 

( 7 

Avg. 

Max. 




Baengnyeong 

849 

0.12 

0.22 

0.98 

606 

0.48 

0.79 

2.79 

0.42 

0.57 

2.65 

Yonsei_Univ (Seoul) 

821 

0.21 

0.30 

1.13 

634 

0.37 

0.66 

2.93 

0.44 

0.35 

1.17 

GwangjuGIST 

1024 

0.16 

0.27 

1.22 

431 

0.40 

0.63 

2.50 

0.30 

0.36 

1.30 

GosanSNU 

967 

0.12 

0.22 

1.14 

488 

0.31 

0.54 

1.95 

0.34 

0.32 

1.47 

Osaka 

1247 

0.11 

0.21 

0.90 

208 

0.22 

0.48 

1.26 

0.14 

0.27 

1.27 

Total 


0.15 

0.24 

1.22 


0.40 

0.65 

2.93 

0.33 

0.40 

1.67 


The calculations were conducted with the assimilated AOD products. The AOD data assimilations were carried out with the GOCI-retrieved AOD products (retrieved via the Yonsei algorithm (Lee et 
al., 2010)) and the CMAQ-simulated AOD products (calculated via Chin et al.’s conversion algorithm (Chin et al., 2002)). 

1 The number of data for non-LRT or LRT events for the entire modeling period; 

2 standard deviation; 

3 AOD increased by LRT events above the average background AOD calculated during non-LRT periods; 

4 the ratios were estimated by ((AODlrt — AOD non _LRr)/AOD non _LRT). Here, AODlrt and AOD non _LRT indicate the average AOD calculated during the LRT and non-LRT events, respectively. 


and Japan was significantly reduced, the stagnant weather 
conditions could seriously deteriorate air quality. Figure 5 
shows a case of the non-LRT events. Noticeable aerosol 
plumes directly transported from China to the Korean penin- 
sula and Japan could not be identified. Some AOD plumes 
were generated in the southern parts of China, and were then 
transported over the northwestern Pacific Ocean. This case 
was classified as non-LRT case since the AOD plumes were 
not directly transported to the Korean peninsula and/or Japan. 

In addition, a high AOD plume (in the second column 
in Fig. 5) was transported through southern Russia (be- 
tween 09:00 and 1L00LST). This plume appears to be thin 
clouds (i.e., cloud optical depth (COD)). This can be con- 
firmed partly via the Terra/MODIS AOD data at 10:30 LST 
on 12 April. The Terra/MODIS AOD data did not show 
any AOD plume over the same region (refer to the third 
and the fourth rows in Fig. 2). Since the GOCI AOD algo- 
rithm has limitations in completely eliminating clouds, thin 
clouds are sometimes considered aerosol plumes. Similar 


cases were also introduced by Al-Saadi et al. (2005) with the 
MODIS AOD data over the United States. Another possibil- 
ity for this plume is that aerosols could be generated through 
Siberian forest fires. In order to check this possibility, the 
global MOZART (Model for OZone And Related chemi- 
cal Tracers) model simulation data was investigated (http: 
//www. acd.ucar.edu/wrf-chem). But no Siberian fire signif- 
icantly affected aerosol plumes on/before 12 April 2011. 

The AOD LRT events were investigated at the five 
AERONET sites (Holben et al., 1998) in the Korean penin- 
sula and in Japan shown in Fig. 1. The AERONET AOD 
values observed at ground sites are regarded as “ground 
truth” because the surface reflectance does not interfere with 
them (Holben et al., 1998; Dubovik et al., 2000; Chin et 
al., 2002, 2004; Park et al., 2011), and due to the high ac- 
curacy (~ 0.01-0.02) of the fully calibrated data (Eck et 
al., 1999). Figure 7 shows the comparison of the GOCI- 
retrieved/assimilated AOD products (red horizontal bars and 
blue solid line in Fig. 7, respectively) with the AERONET 
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(a) CMAQ AOD (b) GOCI AOD (c) Assimilated AOD 

120E 130E 140E 120E 130E 140E 120E 130E 140E 



Fig. 6. As Fig. 3, except that the CMAQ-simulated AOD values 
were calculated via Park et al. (201 l)’s conversion algorithm. 


AOD data (black bars in Fig. 7). The GOCI-retrieved AOD 
values agree well with the AERONET values (e.g., R = 0.84, 
see Table 5). However, both AERONET and GOCI data 
show large time periods with no retrievals, which points 
to the need for models to provide additional information. 
The assimilated AOD product is well correlated with the 
AERONET and GOCI-retrieved AOD values and is much 
more accurate than the a priori CMAQ AOD product when 
compared with the AERONET AOD values. This situation is 
summarized in Table 5, where the errors and biases after the 
assimilation become much smaller than those before the as- 
similation, indicating that the assimilation clearly improves 
the quality of the AOD products. The assimilated product 
provides a continuous estimate of the time series of AOD at 
every location, which the observations alone cannot (the red 
horizontal bars in Fig. 7 only cover limited time periods). 

As discussed previously, the AOD values are related to 
surface PM concentrations. Figure 8 shows an example of 
how to estimate the surface PM concentrations from the 
assimilated AOD values. The daily averaged AOD values 
((AOD)assimiiated) shown in Fig. 7 were converted into PMio 
((PM i o ) estimated ) at five AERONET sites, using following 
proportional relationship (van Donkelaar et al., 2006): 


(PM io) estimated 



x (AOD) assimilated • 0^) 


Twenty-four hour PMio was also measured simultaneously 
at three Korean AERONET locations by the Korean Na- 
tional Institute of Environmental Research (NIER). The ob- 
served 24 h PMio concentrations (black dots) were then plot- 
ted together with the daily averaged PMio concentrations es- 


timated from Eq. (7) (grey bars) in Fig. 8. As shown, there are 
reasonable agreements between the two PMio concentrations 
over two-month periods, except for several time windows 
such as 9-11 April and 13-17 May. This example (Fig. 8) 
is presented here in order to show how AOD is related to 
surface PM concentrations. Although reasonable agreements 
between the observed and estimated PM concentrations were 
shown in Fig. 8, the conversion methods should be further 
developed in order to more accurately convert AOD to PMio 
(or PM2.5) over East Asia, with the consideration of more 
factors, such as height of PBL, RH, and presence of dust and 
BB plumes in the free troposphere (Wang and Christopher, 
2003; Engel-Cox et al., 2004; Liu et al., 2005; van Donke- 
laar et al., 2006; Choi et al., 2009; Schaap et al., 2009; Tsai 
et al., 2011). 

3.2 Estimation of the impacts of long-range transport 

The grey-shaded areas in Fig. 7 represent the periods during 
which the receptors were affected by the PM LRT events. 
The PM LRT events were determined based on the meteoro- 
logical fields generated from MM5 simulations together with 
comprehensive backward trajectory analysis. For example, 
the Baengnyeong site was affected by the small- and large- 
scale LRT events fourteen times over the two-month period. 
Most of the high (and wide) AOD peaks are due to the large- 
scale LRT events from China, and most of the small/sharp 
spikes are usually caused by local pollution under stagnant 
meteorological conditions. The grey-shaded periods in Fig. 7 
usually start and end when the large-scale LRT AOD plumes 
arrive in and depart from the AERONET sites. 

The background average AOD values (calculated exclud- 
ing the AOD values during the LRT events) are represented 
by dashed lines at the five AERONET sites in Fig. 7. The 
average impact of the transboundary PM pollution over the 
entire two-month period (i.e., over “the non-LRT + LRT pe- 
riods”) was calculated using the following formula over the 
two-month period: 

I]/Li{ A( 9Dlrt x XLRT + AOD non _LRT x (1 -Xlrt)} 
AODavg — — » (° 1 


AODinc = A AOD = AODlrt - AOD non _ L RT, (9) 

where AOD avg , AOD non _LRT, and AODlrt represent the av- 
erage AOD values during the entire period and the average 
AOD values calculated during the non-LRT and LRT peri- 
ods, respectively, xlrt is the fraction of the LRT events over 
the two-month period. N is the total number of the grid cells 
being considered (for the calculation at an AERONET site, 
N = 1). After that, the average AOD increase (AODi nc ) dur- 
ing the two-month period was calculated by Eq. (9). The re- 
sults for the analysis for the AERONET sites are summarized 
in Table 4. From this analysis, the average AOD increase of 
117-265 % is found at the five AERONET sites. 
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Fig. 7. Time series analysis of AERONET AOD (black bars) and two AOD products (red horizontal bars indicate the GOCI-retrieved AOD 
products; blue lines represent the assimilated AOD products) over the modeling period of 1 April to 31 May 2011 at five AERONET 
sites: Baengnyeong, YonseiUniv (Seoul), GwangjuGIST, GosanSNU, and Osaka. The averaged background AODs estimated from the 
assimilated AOD values during non-LRT events are presented by blue dashed lines. Also, the periods of LRT events are displayed with 
grey-shaded areas. 


”e 



a. 



Fig. 8. Observed PMp) (black dots) vs. estimated PMp) (grey bars). Daily averaged estimated PMjq was converted from the assimilated AOD 
values at five AERONET sites. 24 h PMp) was measured at three AERONET sites (Yonsei Univ (Seoul), Gwangju GIST, and Gosan SNU). 
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Table 5. Statistical analysis among AERONET, GOCI-retrieved, CMAQ-simulated, and assimilated AODs for the modeling period of 1 April 
to 31 May 2011. 


AOD 


vs. AERONET 


via Chin et al.’s 
conversion algorithm 
(Chin et al., 2002) 

N l 

R 2 

RMSE 3 

MNGE 4 MB 5 

NMB 6 

CMAQ 

946 

0.47 

0.31 

46.52 -0.09 

-14.51 

GOCI 

499 

0.84 

0.19 

42.42 -0.03 

-15.15 

Assimilation 

695 

0.63 

0.28 

35.75 -0.04 

-6.50 

J N: number of data 
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3 RMSE: Root Mean Square Error = E^ (AOD; m — AOD; q) 
4 MNGE: Mean Normalized Gross Error (%) = E^ 


5 MB: Mean Bias = ^ Ef (AOD; M - AOD />0 ) 


m - AOD z,o|\ 

AOD /j0 J 


x 100 


6 MNB: Mean Normalized Bias (%) = L x 100 


The average AOD increases due to the transboundary PM 
pollution over the Korean peninsula were estimated in a sim- 
ilar manner. The Korean peninsula (receptor) was divided 
into three subregions (northern Korea (A), central Korea (B), 
and southern Korea (C)) as shown in Fig. 1 to determine 
subregional differences in the degree of impact from the 
LRT events. The daily variations of the average AOD val- 
ues over the three subregions (A), (B), (C), and the entire 
Korean peninsula are shown in Fig. 9. In panel (a) of Fig. 9, 
CMAQ-simulated AODs were estimated based on Chin et 
al.’s (2002) method, whereas in panel (b), AODs were calcu- 
lated based on Park et al.’s (201 1) algorithm . Panel (b) is pre- 
sented here for the comparison purpose with panel (a). The 
daily and regionally averaged AOD values were presented 
by color codes in Fig. 9. As noticed from Fig. 9, some long- 
range PM transport events affected parts of the Korean penin- 
sula (e.g., 10 and 21 April), while the entire Korean penin- 
sula was also sometimes influenced by the long-range PM 
transport events (e.g., 1 May). The average AOD increases 
(AOD inc ) calculated using Eqs. (8) and (9) in panel (a) of 
Fig. 9 were 128%, 121%, and 112% over the subregions 
(A), (B), and (C), respectively, showing that northern Korea 
experienced the largest impact of the transboundary PM pol- 
lution from China during this period. The overall average es- 
timated AOD increase over the entire Korean peninsula dur- 
ing the two-month period was 121 %. The daily and region- 
ally averaged AOD values in panel (b) of Fig. 9 were in gen- 
eral smaller than those in panel (a), since Park et al.’s (201 1) 
algorithm tends to produce smaller AOD values than Chin 
et al.’s (2002) algorithm, as discussed previously. Again, the 


differences resulting from different algorithms should be fur- 
ther investigated in the future. 

4 Summary and conclusions 

In this study, a new approach was introduced to more 
accurately monitor transboundary FRT events of particu- 
late pollutants over Northeast Asia, combining the CMAQ- 
simulated AOD with the GOCI-retrieved AOD via an OI 
technique (data assimilation). The PM FRT phenomena from 
China (source) to the Korean peninsula and Japan (receptor) 
usually take place by persistent westerlies in Northeast Asia. 
The hourly GOCI-retrieved AOD with high spatial resolu- 
tion of 500 m x 500 m was provided in order to improve the 
accuracy of the CMAQ-simulated AOD through data assim- 
ilation, which produced a two-dimensional set of AOD data 
over the modeling domain. Using the assimilated AOD data, 
the LRT and the non-LRT events were distinguished, and the 
impact of the LRT events were estimated. During the LRT 
events, the average AOD increases of 117-265 % were found 
at the five AERONET sites in Korea and Japan compared to 
the average background AOD value. In addition, the aver- 
age AOD increase of 121 % over the entire Korean peninsula 
was also estimated. This strongly indicates that the notice- 
able amounts of particulate matter emitted from/formed in 
China were transported into the Korean peninsula through 
the LRT processes. 

The analysis framework presented here that utilizes new 
GOCI retrievals of AOD provides a powerful means to bet- 
ter document and estimate the contributions of LRT to PM 
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Fig. 9. Time series plots of daily and regionally averaged AODs, based on AOD products calculated via Chin et al.’s (2002) algorithm (a) and 
Park et al.’s (2011) algorithm (b) over northern Korea (A), central Korea (B), southern Korea (C), and the entire Korean peninsula, shown in 
Fig. 1, for the modeling period of 1 April to 31 May 2011. 


pollution. This has important policy implications. AOD re- 
trievals from geostationary platforms greatly increase the 
amount of observations that go into such estimates and 
greatly reduce the uncertainties that a model adds to the anal- 
ysis. With higher confidence and accuracy in the estimates it 
is easier for policy makers to discuss transboundary PM pol- 
lution over a region like Northeast Asia. Looking forward, 
new observation systems are being designed that will provide 
more comprehensive aerosol information. For example, the 
Korean GEMS and the second version of GOCI (GOCI-2) 
sensors are scheduled to be launched in 2018. GEMS will be- 
come the world’s first geostationary environmental monitor- 
ing sensor that will be able to provide approximately 10 snap- 
shots per day of information on air pollutants over China, 
Korea, Japan, Southeast Asia, and parts of India, Russia, and 
Central Asia, with a grid resolution of 5 km x 5 km. GEMS 
will be equipped with visible and UV channels. Thus, aerosol 
optical properties (e.g., AOD) as well as gas data (e.g., tro- 


pospheric NO 2 , SO 2 , O 3 , CH 2 O, and glyoxal) will also be 
retrieved. 

GOCI-2 will again monitor the ocean colors over a sim- 
ilar domain as GOCI-1, but with 13 visible/NIR channels 
and a smaller grid resolution of 0.25 km x 0.25 km (Ahn, 
2009). Together with smaller grid resolution of GOCI-2, the 
retrieval algorithm of AOPs from the geostationary satellite 
sensors such as GOCI, GOCI-2, and GEMS will be further 
developed in order to overcome the limitations of the cur- 
rent algorithms. For example, as discussed in Sect. 3.1, the 
current retrieval algorithm cannot perfectly screen out thin 
clouds. These limitations and possible errors will/should be 
overcome in the future works. In addition, the AOD dis- 
tributions obtained by applying the methodology discussed 
in this paper with the GEMS and GOCI-2 data will be 
greatly improved. In addition to the data sets that will be 
retrieved from the upcoming satellite sensors (GEMS and 
GOCI-2), the MET and CTM simulations together with the 
emissions of air pollutants in East/Northeast Asia will/should 
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also be improved continuously in order to produce more ac- 
curate AOP distributions. Better AOP distributions will also 
be generated if a better data assimilation technique, such 
as ensemble Kalman filter method, is utilized in the future 
applications. Ultimately, the GEO satellite sensor-derived 
hourly production of more accurate AOPs together with other 
pollutant information over East/Northeast Asia is expected 
to greatly improve our ability to predict and to investigate 
air quality and climate forcing in East/Northeast Asia in the 
future. 

Appendix A 

Over the ocean areas, aerosol retrieval process is hindered 
over sun-glint regions where strong specular reflection oc- 
curs. The sun-glint region is determined by the glint angle 
less than or equal to 40° in both GOCI and MODIS algo- 
rithms, where the glint angle is defined by 

Glint angle = cos -1 (cos 6 cos0° + sinO sin0° cos \j/), 

where 0° is the solar zenith angle, 0 is the viewing zenith an- 
gle, and i// is the relative azimuth angle. Since the sun-glint 
occurs in the specular reflection direction, a larger area is un- 
der the effect of sun-glint for cross-track scanning sensors in 
sun-synchronous orbits, such as MODIS rather than sensors 
in geostationary orbits. 
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